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Mosfet Conclusively, the incremental in ratio has fulfilled the anticipated increment 
Multi gate on the drive current (Ion) as well as reductions of the leakage current (Iorr). 
Short channel effects Threshold voltage (VrH) meanwhile has also achieved the nominal prediction 
that is obliged by the International Technology Roadmap Semiconductor 
(TRS) 2013 for which is at 0.676412.7% V. The Ion, Iorr and Vru obtained 
from the device has evidently met the minimum condition by ITRS 2013 for 

low power Multi-Gate technology. 
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1, INTRODUCTION 

Moore’s law has served as guideline for MOSFET miniaturization for over a decade. Yet, the 
circumstances in the rescaling of the MOSFET is due to the gate controllability decreased over the channel 
caused by high electric field [1, 2]. Besides, shrinking the conventional MOSFET needs improvements to 
bypass the barriers that are due to the fundamental physics [3]. Continuous scaling of bulk-Si MOSFET 
technology have also been massively challenged by the increment in leakage current and also the performance 
variability of transistor [4, 5]. This indicates that the short channel effects (SCEs) have continually been the 
key issues when the rescaling is implemented. That being said, an alternative way to analyse and process multi- 
gate transistor like double-gate FinFET (DG-FinFET) have been done in order to achieve the perfect switch 
for transistor, along with to minimize the SCEs. SCEs occurred due to the charge sharing effect in the shorten 
devices channel, aside from the reason that the channel depletion region and inversion layer channel charge 
control has been contested by both gate and source/drain regions. In addition to that, the extreme increment in 
the subthreshold leakage current has become concerning once it approaches the nanometer regime [6]. 

The electrostatics of the channel has been effectively impacted by the drain potential and, 
consequently, channel has been unable to be controlled by the gate as these are mainly due to the channel 
lengths that are thinned in deeply scaled MOSFET. Consequently, the gate is incapable to fully close the 
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channel in the off-mode operation, which directs to the rise of Iorpr between the drain and the source. 
Meanwhile, the alternative to the existing planar MOSFET, that is the multiple-gate field-effect transistors 
(MGFETs), has showcased a better screening of the drain potential from the channel and this is due to the 
proximity of the supplementary gate(s) to the channel (1.e., higher gate-channel capacitance) [7, 8]. This makes 
the short channel performance metrics, such as subthreshold slope (SS), drain-induced barrier lowering (DIBL), 
and threshold voltage (Vru) roll-off better in MGFETs compared to the orthodox planar MOSFET. These 
metrices improvements would means reduced in degradation in the transistor’s Vry with continuous scaling 
and thus bringing less degradation as well to the leakage current (Iorr) [9]. DIBL effect is evident as the 
shrinkage of the device occurred due to substantial field penetration from drain to source for which the 
aforementioned source and drain depletion regions can infringe the channel without bias due to smaller gaps 
between these junctions in short channel device. The Vrtq roll off in the meantime is caused by reduction in 
Vru due to decrement in gate length as the channel length is unacceptable when the off-state Iorr has become 
too great [10]. Several experiments of designing a 3D FinFET Double-Gate device is done after the initial 
fabrication simulation by revising several physical parameters. DG-FinFET has been contemplated in the 
present work because of its better gate control over the channel which allows an enhanced device 
performance [11]. Figure 1 shows the superior short-channel performance of DG-FinFETs over planar 
MOSFETs with the same channel length [12]. 
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Figure 1. DIBL and subthreshold swing (SS) versus effective channel length for double-gate (DG) and bulk- 
silicon nFETs 


DG-FinFETs emerge superior to planar MOSFETs by overcoming a major source of process 
variation [3]. The FinFET structure also utilizes a fin-shaped body that is perpendicular to the wafer surface to 
carry the current for which it is encased by two of the front and back gates that are developed thinly in order 
to reduce the SCEs. DG-FinFETs do suffer from other process variations. Due to their small dimensions and 
lithographic limitations, DG-FinFETs are subjected to several important physical fluctuations, such as 
variations in gate length, fin-thickness, gate-oxide thickness, and gate underlap [13]. 

Generally, a major amount of work on the projection of leakage current and leakage power prediction 
has been issued by various researchers [11], [14-23]. Bhattacharya and Jha reported that the three-dimensional 
(3D) process simulations and device simulations were carried out so that the fabrication processed and the 
device characteristics of 22nm DG-FinFETs is optimized [3]. Parameters studied has massively enhanced the 
respective Vu, Ion, Iorr values for which is centred on the minimum requirement standardized by International 
Technology Roadmap Semiconductor (TRS) 2013 prediction for low power (LP) multi-gate technology for 
the year 2020. 


2. METHODS AND MATERIALS 
2.1. Virtual Fabrication Process 

In this study, the DG-FinFET fabrication procedure is done through DEVEDIT and ATLAS modules 
from Silvaco International. Both ATLAS and DEVEDIT contrasts in its functionality whereby the ATLAS as 
well as DEVEDIT provides device simulations MOSFET device. The general process flow of a DG-FinFET 
especially in obtaining its electrical properties is contained in the ATLAS module in Silvaco TCAD tool and 
the process simulation methods is summarized as in Figure 2. 
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Figure 2. Basic DG-FinFETs Process Simulation 


Since FinFET have a sizeable outcome towards the device’s Iorr, Vru and Ion, the dimensions of the 
fin has been defined beforehand in order to establish desired effective channel length and also gate width of 
the device [24]. The incremental towards the speed of the transistor switching can be done, should it have a 
high Ion and a low Vry. Based on the constant fin height (Hp) and the fin thickness (Tr), constrain of the 
width can be obtained as: 


W = 2Hpw + Trin (1) 
Leakage performance can also be influenced by the fin shape and its dimension whereby narrower fin 


encourages lower leakage current to be obtain. The values of geometrical parameters are identified as shown 
in Table 1. 


Table 1. Value of the geometrical parameters set for Y-Z plane, X-Y plane and X-Z plane 


Y-Z Plane X-Y Plane X-Z Plane 
Parameters vane Parameters oe Parameters value 
(nm) (nm) (nm) 
Gate Length, Le 30.00 Si0, Width, Wox X-Y 30.0 SiO, Length, Lox X-Z 100.0 
Si0, Length, Lox 100.0 SiO, Thickness, Tox X-Y 30.0 SiO, Width, Wox X-Z 30.0 
SiO. ities 30.00 Metal Gate SiO, Width, Wox 16.0 Metal Gate SiO, Length, 80).0 
X-Y Lox X-Z 
Metal Gate S10, Length, 80.00 Metal Gate SiO, Thickness, 18.1 Metal Gate SiO, Width, 7.98 
Lox Tox x-y Wox x-z 
Metal Gate S102 bs : Aluminium Length, War x. 
Thickness, Ee 2.98 Aluminium Width, War X-Y 10.0 : 10.0 
Aluminium Length, Lar jooo, Ua ucsness eats 15.0 Aluminium Width, Watxz 10.0 
Y 
eine Thickness, 15.00 Polysilicon Thickness, Tp x. 25.0 Polysilicon Leneth, Lawse7 30.0 
AL Y 
ere PHICKUESS: 7.10 Polysilicon Width, Wpm x.y 30.0 Polysilicon Width, Wpu xz 7.0 
DM 
Silicon Thickness, Ts 15.00 Silicon Thickness, Ts x-y 15.0 Silicon Length, Ls x.z 80.0 


Silicon Width, Ws x-y 15.0 Silicon Width, Ws x.z 10.0 
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Six constrained mesh regions were set comprises of silicon as the main substrate, two regions of 
silicon oxide (SiOz), polysilicon, and two regions of Aluminium based bulk bottom electrodes. Meshing 
parameters were then set with its height and width set at 10 nm and 5 nm respectively whereby the maximum 
angle allowed within its triangle is controlled. The maximum slope is bounded at 28° with maximum ratio of 
300. Impurities minimum spacing was applied at 20 nm when X-direction is equivalent to 0. The fabrication 
simulated using geometrical parameters that follows Table | is obtained and analysed from Y-Z plane as in 
Figure 3(A). Therefore, the calibration process was used with standard fluid (Glycerin) which was already 
brought with devices. To validate the data error of the reading from the measurement must less than 0.01. 
Meanwhile, the cross section of DG-FinFET from X-Y plane at elevation = -1 slices the device that showcases 
the Aluminium region and SiO, for which obeys the geometrical parameters in Table 1 as shown in 
Figure 3(B). 

The structure at X-Y plane at elevation = 0 meanwhile exhibits two S102 regions, Silicon and 
Polysilicon for which follows the parameters set as in Table 1 as shown in Figure 3(C). In the meantime, two 
cross sections of X-Z planed DG-FinFET simulated is presented ensuring all regions are viewed in each plane. 
Through geometrical parameters established in Table 1, the structure from elevation = 0.5 in Figure 3(D) shows 
the SiO, region that overlooks Silicon, metal gate SiO2, and Aluminium-based electrode regions contrary to 
the cross section in for when the plane elevation = O, for which veiled the SiO2 region. 
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Figure 3. (A) Cross section of DG-FinFET on Y-Z Plane, (B) Cross section of DG-FinFET on X-Y Plane at 
elevation = -1, (C) Cross section of DG-FinFET on X-Y Plane at elevation = 0, and (D) Cross section of DG- 
FinFET on X-Z Plane at Elevation = 0.5 


2.2. Process Design Fabrication 

For the intent of I-V curves simulation of the devices, the Silvaco TCAD has employed the ensuing 
set of physical models: 1) drift diffusion (DD) model with simplified Boltzmann carrier statistics; 2) Fermi- 
Dirac model as a statistical approach implemented in order to reduce the concentrations of a heavily doped 
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carrier regions. Critical electrical characteristics like threshold voltage, ON-state current, OFF-state current 
and subthreshold swing could be extracted from the Ip vs. Vg characteristic. 

The coefficient of Vps, ALPHA is set at 1.4 1/V alongside the static feedback parameter, GAMMA 
at 0.3 through the application of the ATLAS tool with work function of 4.85 eV. The concentration of p-type 
silicon as the aforementioned main substrate is set at 1x10'® atom/cm? whilst is 1x10?! atom/cm? used for the 
n-type substrate for when Yyin=0, Yuax=0.015, and the characteristic length of laterial in Z-direction is set at 
0.004 at Zyin=0.01, ZmMax=0.0233 and Zmin=0.0767 and Zmax=0.091. Alternatively, the parameters for the 
band-edge alignment is definable as either electron affinity or edge alignment whereby in this design the 
electron affinity is involved. The silicon energy gap at 300K is set at 1.1245 eV followed with 4.05 eV for the 
material affinity at permittivity, e=11.9. The temperature of the bandgap energy equation as well as default 
values applied for each energy bandgap parameters is as following in Table 2 [25]: 


2 
EGALPHA(T, ) 3007 Ty 


E (T,) = E,(0)- "T, +EGBETA EG300 + EGALPHA| 77 opera T, + EGBETA (2) 


Table 2. Default parameters for energy bandgap [26] 


Parameters Value 
Energy gap at 300K, EG300 1.08 eV 
Alpha Coefficient, EGatpua 4.73x107 eV/K 
Beta Coefficient, EGprra 636 K 
Conduction band density, NC300 2.8x10'? cm? 
Valence band density, NV300 2.8x10!? cm? 


Meanwhile, number of times the trap procedure is set to six so that it could be repeated in case of 
discrepancy with the number of initial block iteration set at 45. The block iteration would serve as an iteration 
method alongside Newton when convergence is unattainable through Gummel scheme as the Gummel iteration 
is capable to refine initial convergence guess from Newton method. The Fermi-Dirac model is then 
implemented that adheres the electrons with a semiconductor lattices in such thermal equilibrium at 
temperature T,. The availability electron state is electron-occupied energy € can be obtain through [26]: 


fle) = —— (3) 


, e-f£e 
| + exp[ —— ; 
. kT 
L 


The Fermi level energy in Equation (3) is represented by Er followed with k which denotes the 
Boltzmann constant. Therefore, 8.05 eV is used for when the oxide energy gap is at 300K, with 2.8x10'? cm? 
of conduction band density implemented while the material affinity and its permittivity, ¢ is set at 1 eV and 3.9 
respectively. The polysilicon at permittivity of 11.9 meanwhile requires 1.1245 eV of energy gap at 300K, 4.05 
eV of affinity, and 2.8x10'’ cm? of conduction band density, in addition to the valence band density of 
1.04x10'? cm. Since non-local effects is abandoned in conventional drift-diffusion model of the charge 
transport, the energy balance model is enabled in this model by considering the electron temperatures in the 
Fermi-Dirac statistical method. The parallel electric field dependence model for this design is used as it calls 
on the field-dependent mobility. The maximum magnitude of possible correction is set to 1.0 to damp the 
Gummel iterations. 

In the meantime, the extraction of the drain current against gate voltage (Ip-Vc) and Vu properties 1s 
acquired by applying Newton method when divergence occurred in a solution process whereby decremental of 
multiplication factor is carried out from the initial approximation of the electrode bias step. In this regard, drift- 
diffusion calculations are managed in the. ATLAS using the Newton iteration. 





3. RESULTS AND DISCUSSION 

Based on the 2D structure of the DG-FinFET produced which 1s the cross sections of the device, the 
3D structure is also obtained as shown in Figure 4. Compared to the 2D cross sections obtained, it is observed 
that the 3D structure produced comprises the aforementioned X-Y plane, X-Z plane, and Y-Z plane. The silicon 
region is observed to be concealed away as in Figure 4 due to the fact that it is wrapped and surrounded by the 
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metal gate S102 and polysilicon as presented in Figure 3(C) in forming a notional DG-FinFET structure as this 
has ensured the fin shape minimizes the lIorr, subsequently restricting SCEs from taking its place on the 
structure performance. The variation on gate length, channel doping and S/D doping concentration has been 
implemented to analyze the device performances of Vru, Iorr, Ion and SS and ratio of Ion/Iorr for which a 
curve of Drain Current (A), Ip against Gate Voltage (V), Vas is generated by ATLAS module as in Figure 5. 
Through comparisons made between read value obtained through simulations, against predicted value from 
ITRS 2013, Vrx obtained is slightly higher than predicted value made by the ITRS 2013 which is 0.699 V 
compared to 0.676412.7% V albeit minimal in differences. 
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Figure 4. 3D structure of DG-FinFET 





Nevertheless, the Ion and Iorr both have achieved desired values when assessed with the predicted 
ones whereby the Ion obtained higher than 485 uA/um at 1243.86 uA/um as the Torr attained substantial 
difference in 0.0022 pA/um as compared to 20 pA/um predicted, bringing the Ion/Iorr ratio to be significantly 
higher. Preferably, the device with a greater Ion/Iorr ratio means that the device has met an exceptional value 
due to incremental of Ion versus the decreasing Iorr, that is consistent to the expectation made by the ITRS 
2013 and hence, through Table 3 it is evaluated that the FinFET structure has bring the ratio much superior in 
contrast to the ideal ratio of 24.5 aimed. 


Table 3. Comparisons in value of electrical characteristics based on read value and ITRS 2013 value 


Parameters ITRS Value Read Value 
Vr (V) 0.676412.7% 0.699 
Ton (UA/um) > 485 1243.86 
SS (mV/ dec) 70 ~ 90 71.40 
GIDL - 0 
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Figure 5. Drain Current (A), Ip against Gate Voltage (V), Vas curve of DG-FinFET from ATLAS module 


4. CONCLUSION 

In this study, a 3D double gate FinFET is designed with focus of minimizing the impact of short 
channel effects (SCEs) that consequently contributes to the increment in leakage current (orf) and reduction 
in the drive current (Ion) due to the effects cause by the shrinkage of the structure into the nanometer regime 
that have reduces the channel length. At the end, it is observed that DG FinFET has fulfilled good electrical 
properties with high Ion, and low Iorr based on the electrical characteristics analysed. However, comparatively, 
the VTH does not met the nominal value desired which is differed by 3.3304% higher than, while still within 
the ranges of 0.676+12.7% V predicted by the ITRS 2013. Ideally, the device with substantial difference 
towards the increment of the Ion/Iorr ratio means that the device has met a superlative value that is due to 
increment in Ion and lower Iorr is obtained. In current analysis, the major difference between DG-FinFET with 
conventional planar MOSFET is in terms of leakage current value as it exceeded the minimum value predicted 
which is 20pA as specified by the ITRS 2013. Besides, the device characteristics have met the requirement of 
low power (LP) multigate (MG) technology predicted by the ITRS 2013. 
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